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SUMMARY 


The reduction or alleviation of helicopter vibration will reduce maintenance 
requirements while at the same time increase ride quality and helicopter reliability. 
In forward flight, the helicopter’s fuselage vibration spectrum tends to be dominated 
by multiples of the "N/REV" component. This paper presents a way to use the method 
of adaptive inverse control to identify, in real-time, a controller capable of gen- 
erating N/REV vibration of opposite phase to cancel the uncontrolled N/REV vibra- 
tion component. The control considered in this paper will be that of multicyclic 
feathering of blade pitch. 


INTRODUCTION 


In forward flight, asymmetrical airflow through the rotor disk of the helicopter 
produces a fuselage vibration spectrum which tends to be dominated by multiples of 
the N/REV component; where "N" is the number of blades in the rotor producing 
N cycles of vibration per rotor revolution. The N/REV helicopter vibration 
dynamics may be modeled in the locally linear sense as: 

AZ(j) = T(j)A0(j) (1) 

Here, [Z(j)] is a (p x 1) vector whose elements are the N/REV sine and cosine coef- 
ficients of the Fourier transform of data taken from fuselage accelerometers. The 
vector [ 0 ( j ) ] is an (m x 1) vector whose elements are the selected sine and cosine 
coefficients of the higher harmonic blade pitch. The local linear model postulates 
that within a short time interval the changes in the harmonic coefficients of cyclic 
pitch are linearly related to the changes in the harmonic coefficients of vibration. 
The transfer matrix describing this linear relationship is "T." 

One solution to the vibration control problem would be to identify T, find its 
inverse (for m = p) , and then form the optimal control as: 

A6 (j>opt = -[T(j)] _1 Z(j) (2) 

Z(j) is the total sensed (measured) N/REV vibration. The vector A0(j) is the change 
in multicyclic control necessary to induce a vibration component of opposite phase to 
cancel the sensed vibration. The minus sign thus signifies opposite phase. 

Figure 1 shows this control system. This approach to vibration control is 
limited to transfer matrices which are square, as presented thus far, but can be 
easily generalized for the nonsquare case. 

This paper presents the method of adaptive inverse control as a means to simul- 
taneously identify the inverse of the transfer matrix and to find the vibration feed- 
back gain matrix. This method uses the least mean square (LMS) algorithm of Widrow 


and Hoff (ref. 1) to identify the inverse of the transfer matrix. Thereby, no inver- 
sion is required. Hence, the adaptive inverse control scheme is computationally very 
ef f ic ient . 


REGULATOR DESIGN WITH ADAPTIVE INVERSE CONTROL 


In this paper, the LMS algorithm and the method of adaptive inverse control pre- 
sented by Windrow, McCool, and Medoff (ref. 2) are extended to solve the multiple- 
input/multiple-output helicopter vibration control problem. The necessary input 
command to the controller in this case is a vector whose elements are the uncontrolled 
components of N/REV vibration, changed in sign. In this manner the feedback com- 
mands given to the controller represent the vibration desired to cancel. 

The resultant controller arrangement is shown in figure 2. The controller has 
been placed after the plant to form an adaptation error vector. The commanded 
change in multicyclic pitch at step j is A6(j). Referring to figure 3, it is seen 
that when A0(j) is inputted to the helicopter, a corresponding change in vibration, 

AZ ( j) , results: 

AZ(j) = [T ( j) ] A0 ( j) (3) 

This AZ(j) is subsequently multiplied by C, the controller matrix. If C is the 
inverse of the transfer matrix, T, then [C] * AZ(j) should equal the original com- 
manded change in pitch, A0(j). This is usually not quite equal to the original 
change in pitch because of errors in [C] identification and is thus denoted with a 
'’hat*' symbol, A0(j). 

Adaptation Error = A0(j) - A0(j) 

e ( j) = A0 ( j) - [C ( j ) ] AZ ( j ) 

= A0(j) - [C(j)][T(j)][A0(j)] (4) 

Clearly, this error will be minimized when the controller matrix, C, is the inverse 
of T. This error vector is then used by the LMS algorithm in steepest descent form 
to update the estimate of the inverse transfer matrix, C. This algorithm is an 
iterative identification technique which makes changes in the C (inverse plant) 
matrix elements proportional to the gradient of the expected value of the mean square 
error signal with respect to the current estimate of C, 

C(j + 1) = C(j) + k s V{E[e(j) 2 ]} 

It will be shown that an unbiased estimate of this gradient is given by, 

VE[e 2 (j)] = -2[A0 i (j) - cT ( j) AZ ( j) ] AZ T ( j) (5) 

The steepest descent (LMS-SD) algorithm for calculating each row "i" of the controller 
matrix is then shown to be: 

C*(J + 1) = cT(j) - ZksUO^j) - C^(j)AZ(j)]AZ T (j) (6) 
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LMS ALGORITHM IMPLEMENTATION 


The LMS algorithm (ref. 1) is the iterative estimation algorithm which we will 
use here to identify the inverse of the plant (helicopter) transfer matrix. In this 
section, the LMS algorithm and the method of adaptive inverse control presented by 
Widrow, McCool, and Medoff (ref. 2) will be extended to solve the multiple-input, 
multiple-output helicopter vibration control problem. This extension will be made 
by identifying the inverse transfer matrix in a row-by-row fashion. In this manner, 
the single weight vector originally LMS identified by Widrow and Hoff (ref. 1) now 
becomes a matrix of weights. The LMS algorithm will be adapted to make row-by-row 
corrections to the current row estimates of the inverse matrix. Using the n steepest 
descent” approach, these corrections will be made proportional to the gradient of the 
expected mean square estimation error vector, with respect to the current row esti- 
mate values. In this section we will define the terminology used and explain the 
meaning of the previous statement. 

For each row "i" of the controller matrix at time step "j" we have: 

ejU) = [A^Q) - AZ T (j)C.(j)] (7) 

and the square of this error as: 

e?(j) = [A0 ± (j> - cT(j)AZ(j)][Ae i (j) - AZ T (j)C i (j)] (8) 

Here, A0(j) is the ith element of the control vector fed into the helicopter 
plant (T) . 

The LMS algorithm is an iterative identification technique which makes changes 
in the C (inverse plant) matrix elements proportional to these estimation error 
signals. For each row of the controller matrix, the steepest descent update can be 
written as: 


cT(j + 1) = cT(j) + k s V{E[e?(j)]} (9) 

This equation states that the current row estimate is equal to the previous value of 
the estimate plus a correction term. The correction term is zero when the identifi- 
cation error is zero. 

”E[e|(j)]” denotes the ’’expected value” of the error squared, or the mean 
square error. V{E[e|(j)]} denotes the gradient of the mean square error signal with 
respect to the current estimate of C. That is, 

3{E[e?(j) J> 

V<E[ e |(j)]} ■ 3(;(ltJ) (10) 

An expression of this gradient for the multiple-input, multiple-outout (MIMO) case 
will be presented. Note that if the row error is large, a large change will be made 
to the row estimate of the inverse. If the error is small, a small change will be 
made. The stability constant, k s , governs the stability and rate of convergence of 
this method, as will be shown. 
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The feature which makes the LMS algorithm so efficient is the elaborate manner 
in which the gradient of the mean square error is estimated. This method is now 
discussed in the context of identifying the inverse plant matrix, C. The gradient 

of the mean square error can be approximated for one time step (Ref. 1) as: 

V(E[e?(j)]} = 2e i (j)Vte i (j)] (11) 


Then by noting that the gradient of the error can be expressed as: 

v ei (j) = V[A0 ± (j) - c|(j)AZ(j)] 

= -AZ T (j) 

we may express the gradient of the mean square error as: 

V [ e? ( j ) ] = -2[A6 i (j) - cT(j)AZ(j)]AZ T (j) 


( 12 ) 


(13) 


It is important to realize that in this expression A0(j) is an (n x 1) vector, 
AZ^(j) is a (1 x m) vector, and hence the gradient is not a vector, but an (n x m) 
matrix . 


Although surprisingly simple in form, this approximation to the gradient is 
unbiased. To see this, we note that the square error (for each row i) may be 
expressed as: 

e?(j) = A0?(j) - 2A6 i (j)AZ T (j)C i (j) + c\ ( j) AZ ( j)AZ T ( j)C i ( j) (14) 

The mean square error can then be found by taking the expected value of the error 
squared: 

E[e?(j)] = A0?(j) - 2<D^(Z,0)C i (j) + cJ(j)*(Z,Z)C i (j) (15) 


where the coveriance matrices have been defined as 

<J>T(Z,0) = E[A0 ± (j)AZ T (j) ] (16) 

and 

*(Z,Z) = E[ AZ ( j ) AZ T ( j ) ] (17) 

If we now differentiate the expected value of the error squared with respect to C, 
we find that 

V{E[e?(j)]} = - 2 OT( Z> 0 ) + 2cT (j)4> (Z,Z) (18) 


which is in agreement with the gradient estimate we had before: 
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v [e| (j) ] = -2A0 i (j)AZ T (j) + 2C? ( j) AZ ( j) AZ T (j) 


(19) 


Furthermore, by setting the gradient expression in equation (18) equal to zero, we 
see that each row of the C matrix should converge to the corresponding row of the 
inverse of T in an ordinary least squares sense: 

c ± = 4> -1 (z,z)$T(z,e) 


= [AZ(j)AZ T (j)]- 1 A0 i (j)AZ T (j) 


( 20 ) 


as this is the same solution as would be obtained for the ordinary least square error 
solution. 

Thus by using 


V[e*(j)] = -2[A9 i (j) - cT(j)AZ(j)]AZ T (j) (21) 

as the gradient estimate of the mean square error, the LMS-SD algorithm for each row 
"i" of the controller matrix is: 

C*(j + 1) = C^(j) - 2k s [A6.(j) - cT( j)AZ (j) ] AZ T (j) (22) 


Stability and rate of convergence are determined by the selection of the stabil- 
ity constant, k g . Considerations governing the selection of k s will now be 
discussed. 


SELECTION OF STABILITY CONSTANT k g FOR ALGORITHM STABILITY AND CONVERGENCE 


It has been shown thus far that we can find an unbiased estimate of the mean 
square error gradient for steepest descent utilization. The convergence properties 
of the steepest descent technique will now be discussed. Specifically, there is a 
need to know when the method can be guaranteed to converge to the true inverse of 
the plant, and how long it will take to do so. Our analysis will follow that of 
Widrow and Hoff (ref. 1) and will be extended to incorporate the helicopter (plant) 
inverse identification problem. 

For each row (i) of the C (inverse plant) matrix the following difference equa- 
tion was written for the method of steepest descent: 


C^(j + 1) = C^(j) - 2k s [A0 1 (j) - (j ) AZ (j) ] AZ T ( j ) 

= C*(j) - 2k s A6 i (j)AZ T (j) + 2k s cT(j)AZ(j)AZ T (j) (23) 


Taking the expected value of both sides of this equation we get: 
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E[C?(j + 1)] = [I + 2k s <KZ,Z)]E[cT(j)] - 2 k s $T(Z, 0 ) 


J 

= [I + 2k s $(Z,Z)]J +1 C(0) - 2k g £ [I + 2k s $(Z,Z)]<i>?(Z,e) (24) 

i=o 

This equation can be put into modal form by writing the signal covariance matrix as 

*(Z,Z) = Q -1 AQ (25) 

where Q is the matrix composed of the eigenvectors of $(Z,Z) and A is a diagonal 
matrix whose elements are the eigenvalues of $(Z,Z). Substituting this into the 
above equation yields: 

j 

E[C?(j + 1)] = [1 + 2k s Q _1 AQ]j +1 C i (0) - 2kg ^ (I + 2k s Q _1 AQ) (Z , 6) 

i=o 

(26) 

j 

E[cJ(j + 1)] = Q _1 [I + 2k s A]j +1 QC^(0) - 2k s Q“ 1 E (I + 2k g A) *<}*£ (Z , 6) 

i= o 


From this we note that as long as the elements of 

[I + 2k s A] 


are all less than one, then the term 

Q -1 [ I + 2k s A] j +1 Qc|(0) 

will converge to zero, and the second term will converge to 

j oo 


limit - 2k s Q 1 Y 
j-x» i= o 


[I + 2k s A] i Q^(Z,0) = -2kgQ _1 


E (1 + 2k s A p ) 1 Q<t|(Z,0) 

i=o 


2k s Q [\ - (1 + 2kgX p )j Q'J’i (Z > 0 ) 
= " 2 k s Q ~ 1 [2^] Q $T i( z * e > 


= Q -1 A x Q<i>T(Z,e) 

= (Z , Z) (z , 0) 
= C LMS 


(27) 
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Thus, for convergence to the correct solution by the method of steepest descent, 
we see that 


1 + 


2k sVaxl < 


(28) 


T~ < k s < 0 (29) 

A max 

Within these limits, the value of k s is selected with regard to the desired 
application. Values of k s close to zero will be less apt to be affected by random 
noise variations on the plant at the expense of slow adaptation. Values of k s near 
-1/A max will adapt rapidly, but will be more prone to tracking random noise after 
"convergence" has been achieved, and will tend to oscillate about the correct solu- 
tion. The right value of k s is the one which converges at a sufficiently rapid 
rate, yet does not produce unacceptable steady state convergence errors. It is pos- 
sible to make the value of k s depend exponentially on the time step, but this will 
not be discussed here. 

The rate of convergence may be thought of in terms of learning curve time con- 
stants, t. Letting rp denote the geometric ratio of the pth mode, we can say 
that 


r p = (1 + 2k s X p ) 


(30) 


Then, 


or 



2 M P 


(31) 


(32) 


If all eigenvalues of $(Z,Z) were the same, a single exponential learning curve 
could be defined. In the more general case, however, the eigenvalues will not 
be equal. Then the learning curve will be a function of all of the eigenvalues 
corresponding to the various normal modes. Thus, we can expect the fast modes to 
cause rapid initial learning, whereas the learning errors associated with the slower 
modes will take longer to die out and govern final convergence. 


STARTING THE ADAPTIVE INVERSE VIBRATION CONTROL ALGORITHM: THE TRAINING PHASE 


We have thus far established a method for iteratively correcting the estimate 
of the helicopter transfer matrix inverse. Yet, nothing has been said about initial 
conditions. At the start of the control routine, whether on the ground or in the 
air, a good starting value for the inverse of the helicopter transfer matrix may not 
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be known. Hence, a learning phase, termed the "training phase," is required at the 
start of the algorithm. 

During the training phase, the blade pitch is given small perturbations and the 
corresponding changes in vibration are measured. These measurements are then used 
by the LMS-SD algorithm to correct an arbitrary initial value of the inverse plant 
matrix, C(0). Any value of C(0) may be used — for example, all elements equal to 
zero. No vibration control commands are generated during the training phase. In 
this manner, large transients in control are avoided. Hopefully, after a number of 
iterations, the inverse plant matrix will be accurate enough to use in the actual 
control phase. 

The control phase begins when the sum of squares of the adaptation errors, 

A0 - A0, becomes small enough. No further training is then required. The LMS-SD 
algorithm will update the estimate of the plant inverse transfer matrix to keep up 
with changes in the plant (helicopter) or changes in the operating environment. 


SIMULATION OF A THREE-INPUT, THREE-OUTPUT PLANT 


The LMS-SD algorithm presented in this paper was simulated for a three-state 
plant. An arbitrary three-by-three transfer matrix, T, was chosen to represent the 
helicopter (fig. 4). The helicopter-vibration dynamics were then modeled in the 
global sense by the following equations: 


"z 1 (jf 


' T n 
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T 1 

13 
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= 
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The LMS-SD update equations which identify the inverse plant transfer matrix (con- 
troller matrix) were then: 

c?(j + 1) = C? (J ) - 2k s [A9 1 (j) - cT(j)AZ(j)]AZ T (j) 


c 2 (j + 2 ) = c 2 (j) 


2k s [A0 2 (j) - C 2 (j)AZ(j) ]AZ T (j) 


(34) 


C|(j + 3) = C?(j) 


2k s [A0 3 (j) - cT(j)AZ(j)]AZ T (j) 


And the vibration control commands were then generated by: 


■01 u 

+ 

If 


■01 (jf 


C 11 

Cl2 

^13 


"AZ^jf 

0 2 (j 

+ 

1) 

= 

e 2 (j) 

- 

C 21 

C 2 2 

C 23 

* 

AZ 2 (j) 

_0 3 a 

+ 

1). 


_0 3 ( j )_ 


- C 3 1 

^3 2 

^3 3- 


_AZ 3 (j)_ 


( 35 ) 
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The simulation included both the training and the control phases. The effects 
of identification errors and changes in the operating conditions were also simulated. 
The overall program is diagrammed in figure 6. 

From the preliminary simulation results, the LMS-SD algorithm appeared to pro- 
vide a very robust control system. Figures 4 and 5 illustrate some of the interest- 
ing results obtained. Figure 4 illustrates the typical convergence pattern for k s 
chosen near its midrange stability. Here the sum of the adaptation errors, A0 - A0, 
have been plotted versus the identification iteration cycle count. From this plot 
the normal modes of the system can be seen. Initial rapid convergence seemed to be 
produced by the fast modes, whereas final convergence appeared to be governed by the 
slower modes. 

Figure 5 shows the behavior of LMS-SD convergence and vibration levels through 
the training and control phases. To simulate a sudden change in the operating condi- 
tions, the plant simulation values were suddenly changed by 20% and the uncontrolled 
vibration level was changed by 40%. The effect of these changes is shown in figure 5. 
The controller was able to respond in a manner to alleviate the vibration and a 
rather surprising feature of this robust control system was discovered. 

The controller values did not have to converge exactly, or even very closely, to 
cancel the vibration. For example, figure 5 shows that after a step disturbance in 
the uncontrolled vibration level, the controller may retain steady-state identifica- 
tion errors, while rapidly alleviating the vibration. The explanation of this robust- 
ness rests in part upon the identification and control phases being iterative pro- 
cesses. As such, the proper change in blade-pitch control can be generated by a 
number of linear combinations of control changes over several steps. It turned out 
(in simulation) that when the controller matrix was sufficiently close to the true 
inverse, a linear feedback combination capable of alleviating most of the vibration 
was generated. When the vibration level approached zero, the pitch vector was no 
longer altered, thus preventing complete identification of the inverse. This was as 
it should be, since if there is no vibration, there is no need to change the pitch- 
control vector. The controller matrix can then no longer be altered since the update 
correction term (eq. (22)) goes to zero for no uncontrolled vibration. In actual 
implementation, however, random noise and perturbations would most likely allow com- 
plete convergence to the true inverse. The significant point though is that complete 
convergence is not required for good vibration control. 


CONCLUSIONS 


The method of adaptive inverse control was presented in the context of helicopter 
vibration control. The theory used to derive the MIMO LMS-SD algorithm was given. 
Results of computer simulations of the inverse adaptation process were also presented 
and shown to indicate the validity of the results predicted by the theory. The con- 
clusions verified by simulation were as follows: 

1. The method of adaptive inverse vibration control can provide a very robust 
control system which does not require a perfect knowledge of the inverse plant matrix 
to alleviate vibration. 

2. The "training phase" and the LMS-SD algorithm allows the method of adaptive 
inverse control to be implemented without any a priori knowledge of the plant. No 
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dynamic models are required, as the training phase can quickly determine proper con- 
troller initial conditions. 

3. The LMS-SD algorithm is capable of adapting the controller to track changes 
in the flight conditions or changes in the uncontrolled vibration level. The con- 
troller converges quickly for moderate disturbance levels, while taking longer for 
larger (step) disturbances. 

4. The learning curve of the controller during the training phase was shown to 
be quantitatively close to that predicted by the learning curves of the normal modes. 

5. The method of adaptive inverse control should have very fast on-line capa- 
bility, as implementation of the LMS-SD algorithm requires so few operations. 
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Figure 3.- Formation of the adaptation error used by the LMS-SD algorithm to 
correct the estimate of the inverse plant matrix, C. 


13 


k $ = -0.30 


1.5 


T = 


1 0 
3 1 

1 2 


SUM SQUARE 
IDENTIFICATION 
ERRORS IN C 



Figure 4.- The "learning curve" of the adaptation process, or a plot of the sum 
square error in the controller (inverse) identification versus the number of 
iterations from initial conditions; k s chosen near the middle of its stability 
range . 
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Figure 5.- Controller convergence properties during the training and control phases, 
along with the uncontrolled vibration. The effect of a 20% change in the plant 
transfer matrix, T, and a 40% change in the uncontrolled vibration. 
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Figure 6.- LMS steepest descent program. 
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